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ABSTRACT 

In this paper we define an observationally robust, multi-parameter space for the classification of nearby 
and distant galaxies. The parameters include luminosity, color, and the image-structure parameters: 
size, image concentration, asymmetry, and surface brightness. Based on an initial calibration of this 
parameter space using the "normal" Hubble-types surveyed by Frei et al. (1996), we find that only 
a subset of the parameters provide useful classification boundaries for this sample. Interestingly, this 
subset does not include distance-dependent scale parameters, such as size or luminosity. The essential 
ingredient is the combination of a spectral index (e.g., color) with parameters of image structure and 
scale: concentration, asymmetry, and surface-brightness. We refer to the image structure parameters 
(concentration and asymmetry) as indices of "form." We define a preliminary classification based on 
spectral index, form, and surface-brightness (a scale) that successfully separates normal galaxies into 
three classes. We intentionally identify these classes with the familiar labels of Early, Intermediate, and 
Late. This classification, or others based on the above four parameters can be used reliably to define 
comparable samples over a broad range in redshift. The size and luminosity distribution of such samples 
will not be biased by this selection process except through astrophysical correlations between spectral 
index, form, and surface-brightness. 

Subject headings: galaxies: classification — galaxies: morphology — galaxies: colors 



1. INTRODUCTION 

It is now well-established that a large fraction of galaxies 
discovered at intermediate and high redshift have unusual 
morphologies, and thus cannot be classified in terms of 
the nominal Hubble-Sandage system (Driver et al. 1995, 
1998; Abraham et al. 1996a, 1996b). The Hubble clas- 
sification scheme is also difficult to apply to many local 
galaxies, dubbed 'peculiar,' or any galaxies imaged at low 
signal-to-noise (S/N), or apparently small size (relative 
to the point-spread function). The Hubble-Sandage clas- 
sification system was predicated on the study of nearby, 
"nor mal" galaxies - luminous and relatively quiescent ob 
iectsf Bandage 1961| , |5andage 



Tamman 1987 
Bcdkc 1993 ). While the classification system 
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de Vaucouleurs et al. (1976) makes an attempt to push the 
framework to "later" types, it still suffers from the above 
shortcomings. Fundamentally, these traditional classifica- 
tion schemes are based on the concept of 'pigeon-holing' 
galaxies based on a reference-set, or archetypes. These 
archetypes are selected from samples in the local universe, 
and are preferentially axisymmetric systems. Since our lo- 
cal census is undoubtedly incomplete, and, since galaxies 
evolve, such reference sets by their very definition are in- 
complete. Thus it is not surprising that these systems are 
of marginal utility in the study of dwarf galaxies, interact- 
ing galaxies, or galaxies at high redshift. 

An alternative classification scheme could be based on 
quantitative indices, the inter-relation of which is not pre- 
determined by a finite reference set. This would permit 



galaxies to be classified, for example, in different stages 
of their evolution; albeit the classification would be differ- 
ent but the basis set of indices would be the same. The 
goal of this paper is to define such a set of indices that 
can be used as quantitative, objective classifiers of galax- 
ies (i) over a wide range in redshift, and (ii) for wide range 
of galaxy types. In particular, we desire classifiers that 
are well suited to typing both "normal" galaxies and the 
compact galaxies that are the focus of a companion study 
(Jangren et al. 2000; hereafter, paper II). We anticipate 
that such a classification scheme is both necessary and 
enabling for the exploration of the physical mechanisms 
driving galaxy evolution (Bershady 1999). 

What are the desired characteristics of classification pa- 
rameters? They should be physically interesting (closely 
related to underlying physical properties of galaxies), 
model- independent, and measurable for all galaxy types. 
It also should be possible to accurately determine the pa- 
rameters chosen for a wide range of image resolution and 
signal-to-noise ratios. 

From Hubble's classification a posteriori we have learned 
that a strong correlation exists between galaxy spectral 
type and apparent morphological features - at least for the 
galaxy types which fit well within his scheme. This corre- 
lation - noted by Hubble as early as 1936 (Hubble 1936) 
- can loosely be termed a 'color-morphology' relation, al- 
though the correlation is not necessarily limited to broad- 
band color. This is a triumph of Hubble's classification 
explicitly because it is not part of the classification. Fur- 
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thermore, the correlation yields clues about the physical 
connection of the present matter distribution and the star- 
formation histories in galaxies. But while morphology (or 
'form') and spectral type are correlated, there is also signif- 
icant dispersion in this correlation. Some of the more no- 
table deviations from the nominal color-morphology rela- 
tion are found in the plethora of forms for spectrally 'late' 
type galaxies, the presence of 'E+A' galaxies (Dressier & 
Gunn 1993), and the compact, luminous, blue, emission- 
line galaxies studied in paper II (as we shall show). This 
points to the importance of form and spectral type as key, 
yet independent axes of a revised classification system. 

However, the only example of such a revised classifica- 
tion system is that of Morgan (1958, 1959), where central 
light concentration is used as the primary classification 
parameter. Morgan was motivated by the fact that (i) a 
salient criterion used in classifying galaxies in the Hubble- 
Sandage system is the degree of central concentration of 
light; (ii) there wa s a significant dispersi on in spectral type 
and Hubble type (Humason et al. 1956); and (iii) spectral 
type appeared to correlate more strongly with light con- 
centration. In this way, Morgan hoped to wed the classifi- 
cation of stellar populations to the classification of galax- 
ies. Nonetheless, he was compelled to introduce a sec- 
ondary parameter, i.e. the 'Form Family' because there 
was still a dispersion of morphological forms within each 
of his spectral types. Today, one should be able to im- 
prove upon Morgan's scheme by introducing quantitative 
measures of image concentration and other indices of form, 
and by independently assessing the spectral type via colors 
or spectra. 

A number of subsequent attempts have been made to 
construct quantitative classification system that could re- 
place or modify the current Hubble scheme. Yet these 
schemes are ge nerally based purely either on photomet- 
ric fo rm (e .g. Elmcgrccn fc Elmcgrccn 1982|; pkamura 
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Watanabe et al. 1985: Poi et al. 1993: Abra- 



1994|; Pdcwahn 1995|; Han 1995|) or spectra 
Bershady 19951; IConnolly et al. 19951; |Zaritskv 



19 951; Folkes et al. 1996| ; |Bromlcy et al. 1998| ; Ho, 
. 1999[ ) ■ I 11 essence, they have relied implicitly on 



an assumed correlation between galaxy spectral type and 
apparent morphology. Related attempts have been made 
to use artificial neural networks to reproduce the Hubble 
scheme in an objective way (e . g. |Burda &: Fcitzinger 1992 



et al. 1993J ; |Naim et al. 1995| ; pdcwahn 1995| ; |Odc"l 
et al. 1996| ). Yet these go no further in diffcrentiat- 
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ing between spectral type and form. Only in Whitmore's 
(1984) scheme are spectral and structural parameters com- 
bined, i.e., B — H color, size, and bulge-to-total ratio are 
used to define two principal classification axes of scale and 
form. But again, the correlation(s) between galaxy spec- 
tral type, scale and form are not explicit. 

Here we attempt to expand on Morgan's program by 
fully quantifying the classification of form via image con- 
centration and several other structural parameters, and 
explicitly using color as an indicator of spectral type. In 
this study we choose to use only a single color (B — V), but 
we anticipate that a more desirable, future development 
would be to include broad-wavelength coverage, multi- 
color data and spectroscopic line-indices. Spectroscopic 
line-indices would be required, for example, to identify 



E+A galaxies. While such galaxies are not the focus of 
the present work, a comprehensive classification scheme 
should be able to isolate these systems and determine the 
range of their morphology (cf. Dressier & Gunn 1992, 
Couch et al. 1994, and Wirth et al. 1994). Nonetheless, 
broad-band colors are a cost-effective way to characterize 
the spectral continuum (cf. Bershady, 1995, and Connolly 
et al. 1995). Of more direct relevance to the study at 
hand, a future elaboration of including U — V and V — K 
would enhance the ability to distinguish between spectral 
types particularly for galaxies with extremely blue, optical 
colors (e.g. Aaronson 1978, Bershady 1995). 

We have also chosen to quantify form and scale via 
non-parametric measures, such as luminosity, half-light 
size and surface-brightness, asymmetry, and image con- 
centration. An alternative, model-dependent approach 
is to decompose a galaxy's light profile into a disk and 
bulge. The traditional one-dimensional decompositions 
are fraught with technical problems such that decompo- 
sitions can on ly be achie ved reliably for about half of all 
disk galaxies (Kent 1985). The newer two -dimensional de - 
composition techniques are superior (e.g. de Jong 1996b), 
and have been shown to successf ully reproduce obse rved 
light profiles for faint galaxies (e.g. Simard et al. 1999| ). In- 
deed, one can argue that two-dimensional model fitting to 
imaging data is optimum in terms of using the available in- 
formation, and for minimizing random error. At high S/N 
and high angular resolution, however, even the most "nor- 
mal" galaxies exhibit peculiarities (as discussed in more 
detail in §3.3.2) such that simple bulge-plus-disk models 
cannot reproduce these frequently observed peculiarities in 
light distributions with high fidelity. The situation wors- 
ens for "peculiar" galaxies. For this reason we have some 
concerns about the uniqueness of the observationally de- 
rived model parameters, and hence their interpretation. 
We anticipate future developments which use the models 
and non-parametric measurements in a hybrid scheme op- 
timal for characterizing galaxy light distributions both in 
terms of random and systematic errors. 

It is worth noting again that bright galaxy samples are 
notorious for missing or under-representing certain galaxy 
types - particularly dwarfs and low-surface-brightness 
galaxies. The samples used here are no exception. While 
this was one of our complaints about the classical Hub- 
ble scheme, there are two key differences with our ap- 
proach: (i) the classification parameters we develop are 
objective; and (ii) these parameters do not assume the 
presence of basic axi-symmetry, disk-plus-bulge structure, 
or spiral patterns which underly the Hubble scheme. As 
we will show, the galaxies examined here are sufficiently 
diverse to establish the parameter space for a comprehen- 
sive classification scheme, although not the comprehensive 
classification itself. By developing an initial classification 
of these galaxies, however, we intend to use it as a foil 
against which we can begin to compare the classification 
of more distant samples: How are the classifications differ- 
ent? Do the nearby and distant samples occupy the same 
regions of parameter space? If not, do the differences rep- 
resent continuous extensions of these parameters, or are 
they physically disjoint? These are the types of questions 
one can address given the limitations of current local sam- 
ples. Note that we must stop short of identifying differ- 
ences as "new," epoch-specific classes of galaxies. Without 
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a complete census of both the nearby and distant universe, 
it is not possible to establish whether there are different 
"classes" of galaxies at different redshifts; apparent dif- 
ferences could simply be artifacts of the presently limited 
samples. With such a complete census, in the future we 
may hope to address deeper issue of how the comoving 
space-densities of different classes evolve. 

Towards the goal of establishing a comprehensive classi- 
fication scheme of utility to distant galaxy studies, in this 
paper we assemble a robust set of non-parametric, pho- 
tometric and structural properties for a range of nearby, 
lumimous galaxies. We define a multivariate, photometric 
parameter space that forms an initial classification scheme 
for these galaxies. This classification can be used reli- 
ably to identify comparable samples in other surveys and 
at higher redshift. In the accompanying paper (II) we 
measure these properties for compact, luminous emission- 
line galaxies at intermediate redshift, compare them to the 
"normal," nearby galaxies studied here, and demonstrate 
that our classification parameter space distinguishes be- 
tween these two samples. We discuss the implications for 
the evolution of this intermediate-redshift sample therein. 
In future papers in this series we intend to extend our 
analysis (a) to more representative samples of the local 
volume that include dwarf and emission-line galaxie s (e.g 



the U niver sity of Michigan Objective Prism Survey (Salzei 



et al. 1989); (b) to more comprehensive samples of distant 
galaxies, e.g magnitude-limited samples from the Hubble 
Deep Field; and (c) to studies of the morphological evo- 
lution of these distant samples. The classification scheme 
which we propose here is intended as a framework for these 
future studies. 

The data sets are presented in §2; the analysis is de- 
scribed in §3. The results are presented in §4, and sum- 
marized in §5. Throughout this paper we adopt Hq = 50 
km s" 1 Mpc~\ q = 0.1, A = 0. 

2. NEARBY GALAXY SAMPLES 

As a primary reference sample, 101 of the 113 local 
Hubble-type galaxies from the catalog of Frei et al. (1996) 
were analyzed. This sample will define what we mean by 
"normal" galaxies in this paper. This catalog is the only 
digital, multi-band, sample publicly available that is rea- 
sonably comprehensive; it consists of ground-based CCD 
images of bright galaxies, all apparently large (most have 
diameters of 4' - 6') and well resolved. As a result, the sam- 
ple contains mostly luminous and physically large galaxies: 
out of the 101 objects we used in our analysis, only seven 
have L < 0.1 L*. We excluded 12 objects whose apparent 
sizes were larger than the CCD field of view (thus their 
image structure parameters could not be well estimated). 
Two of the excluded objects are early-type galaxies (E- 
S0), seven are intermediate (Sa-Sb), and three are late- 
type (Sc-Irr). 1 The majority of the remaining sample are 
spirals and SO galaxies. Frei et al. have removed fore- 
ground stars from the images of the nearby galaxies, in 
a few cases leaving visible "scars;" except in the case of 
NGC 5792, these residuals did not cause noticeable prob- 
lems when determining the structural parameters (§3.3). 

In several instances in the present analysis we reference 
the sample of Kent (1984, 1985), which is composed of 53 



nearby, luminous and physically large galaxies similar to 
the Frei et al. sample. We find Kent's sample useful for 
comparison of both photometric and structural parame- 
ters. We also reference the sample of 196 normal (non- 
active) Markarian galaxies studied by Huchra (1977a). 
Relevant characteristics of the above three samples are 
summarized in Table 1 , including an enumeration of the ef- 
fective filter systems used in each study. Further details on 
these photometric systems are found in the studies listed 
in the Table and references therein. 

2.1. Comparison of reference samples to emission-line 
galaxy samples 

Both the Frei et al. and Kent samples are under- 
representative of dwarf galaxies, and contain neither HII 
galaxies nor low surface-brightness galaxies. The latter ob- 
jects have been shown to make up a significant fraction of 
the local galaxy population (de Jong 1995, 1996a). Clearly 
our reference samples do not constitute a representative 
template of the local population. Here we estimate where 
these samples may be particularly un-representative with 
an eye towards the study of faint galaxy samples in future 
papers. In Figures 1 and 2 we compare the Frei et al. sam- 
ples photometric properties of color and luminosity to (i) 
the normal Markarian galaxies (Huchra 1977a), (ii) dwarf 
spheroidals (as described in the following section), and (iii) 
the intermediate redshift samples presented in paper II. 

Since the Markarian galaxies were selected from objec- 
tive prism plates based on their strong UV continua, the 
sample is biased toward bluer colors than the Frei et al. 
galaxies and is thus likely more representative of star- 
forming galaxies. Huchra's sample contains fainter galax- 
ies that extend the magnitude range down to Mb ~ — 14 
and the color-color locus blue-ward of B — V = 0.4. 

The intermediate-redshift galaxies, also selected in part 
due to their blue color (see paper II), have blue lumi- 
nosities comparable to the brighter half of the Frei et al. 
sample, but with bluer colors. This places most of them 
in a distinct region of the color-luminosity plot from the 
Frei et al. sample. In contrast, the distribution of the 
Markarian galaxies extends into the region occupied by the 
intermediate-redshift objects. In the color-color diagram, 
again the intermediate-redshift galaxies largely overlap 
with the Markarian sample in the region corresponding 
to extreme blue colors not occupied by the Frei et al. or 
Kent samples. 

In short, the Frei et al. sample is spectro- 
photometrically disjoint from extreme samples of blue, 
star-forming galaxies at intermediate-redshift (e.g., paper 
II), even though both contain intrinsically luminous and 
moderate-to-high surface-brightness systems. Yet clearly 
there ar local examples (e.g., from Markarian) which are 
as blue and luminous as these intermediate-redshift, star- 
forming galaxies. These sources are simply missing from 
the Frei et al. sample. The comparison of the global prop- 
erties of the intermediate-redshift, compact, star- forming 
galaxies in paper II to those of local galaxies from Frei 
et al. (here) is then an initial step in mapping the range 
of galaxy types at any redshift. Further investigation 
of the nature and evolution of these types of extreme, 
star-forming systems will be greatly facilitated by future 



lr The excluded objects are: NGC 2403, 2683, 3031, 3079, 3351, 3623, 4406, 4472, 4594, 4826, 5746, and 6503. 
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work quantifying the image structure of local counterparts 
B-V< 0.4 and M B < -19. 

2.2. Comparison of reference samples to dwarf 
spheroidals 

We have made some attempt, where possible, to ac- 
cess the photometric and structural properties of other key 
dwarf populations. We schematically indicate the locus of 
dwar f ellipticals/spheroidals in Figures 1 and 2 us ing data 
from jCaldwcll 1983 , Bingelli & Cameron 1991, and Bingellj 
& Jerjen 199S . The dwarf spheroidals occupy a virtually 
unpopulated region of the color-luminosity diagram at rel- 
atively red colors and low luminosity. The absence of such 
objects from most surveys is attributed typically to a selec- 
tion bias since these sources are at low surface-brightness. 
It is interesting to note that in the color-color diagram the 
dwarf spheroidals occupy a region over-lapping with the 
early-to-intermediate type spirals. Hence the integrated 
broad-band light of these systems are unusual compared 
to our reference samples only with respect to their lumi- 
nosity. We refer to the dwarf spheroidal properties exten- 
sively in future papers where we also explore their image 
structural properties. 

3. ANALYSIS 

As noted in the Introduction, many galaxies are suf- 
ficiently unusual that they cannot be classified in terms 
of the normal Hubble scheme. This becomes increasingly 
true at intermediate redshifts. The compact, luminous 
emission-line galaxies in paper II are such an example. 
This is not due to poor spatial resolution, but to truly un- 
usual morphological properties, e.g., off-centered nuclei, 
tails, asymmetric envelopes, etc. To compare such objects 
morphologically to "normal" galaxies, we define here six 
fundamental parameters of galaxy type that are quantita- 
tive, can be reliably determined over a range in redshift, 
and are physically meaningful. 

Two of these parameters are photometric, derived 
from existing ground-based imaging and estimated k- 
corrections: rest-frame color B — V, and absolute blue lu- 
minosity Mb- Two are image structure parameters, de- 
rived from multi-aperture photometric analysis presented 
below: physical half-light radius R e , and image concentra- 
tion C. One is a combined photometric-structural param- 
eter: average rest- frame surface brightness SB e within R e . 
Of the three parameters luminosity, half-light radius, and 
surface brightness, any one can obviously be derived from 
the other two. (We consider all three since in any given 
range of, e.g., luminosity, there is significant dispersion in 
both SB e and R e .) The sixth parameter, a 180°-rotational 
asymmetry index (A) , utilizes the multi-aperture photom- 
etry indirectly through definition of the extraction radius 
for rotation; we refer to A as a structural parameter. Table 
2 contains all individual measurements for the Frei et al. 
sources. Luminosities and all image-structure parameters 
are measured in the rest-frame B band. 

3.1. Photometric parameters: restframe color and 
luminosity 

While the Frei et al. (1995) data set contains Bj images 
for 75% of the sample and g band images for the remaining 
objects, there is no blue bandpass in which observations 



are available for all galaxies (see Table 1). We have made 
a comparison of the apparent (uncorrected) B magnitudes 
l isted in the Third Reference Catalogue of Bright Galaxies 
(de Vaucouleurs et al. 1991, RC3) to those derived from 



our photometry of the Frei et al. images, appropriately 
transformed to the -B-band using the tabulated correc- 
tions of Frei and Gunn (1994). This comparison shows 
that while the two magnitude estimates do not differ in 
the mean, there is a 0.25 mag (rms) scatter. To avoid the 
uncertainty associated with SED-dependent color trans- 
formations (see also §A.3) we use the RC3 uncorrected B 
magnitudes and B — V colors instead of the values from 
our own photometry. 

We apply fc-corrections and corrections for galactic ex- 
tinction to the B — V colors and apparent B band magni- 
tudes of the nearby galaxy sample in the manner described 
in RC3. The heliocentric velocities V) lc of the galaxies are 
small (no greater than 3000 km s _1 for any object); the 
average velocity is (vh c ) ~ 1000 km s _1 . Hence the asso- 
ciated fc-corrections are < 0.05 mag. We use the dista nces 
given in the Nearby Galaxies Catalogue ( Tully 1988| ), re- 



calculated to f/o=50 km s Mpc , to derive absolute 
magnitudes Mb from the corrected apparent magnitudes. 
Note we do not correct for internal extinction since the 
suitability and procedure for applying such corrections 
may be ill-defined for higher-redshift galaxies. 

3.2. Multi- aperture photometry 

To characterize the light distributions of the galaxies, we 
performed multi-aperture photometry on all images. The 
apertures are centered at the intensity-weighted centroid 
of each object. Since much of the profile shape information 
is contained in the central parts of the image, logarithmi- 
cally spaced apertures are used. For the photometry of 
the Hubble Space Telescope (HST) images in paper II, 
the smallest aperture corresponds to 0."05, the largest to 
15"; for the nearby galaxy photometry here, the aperture 
radii are scaled to correspond to similar linear sizes. The 
apertures are circular to accommodate the irregular mor- 
phology of the intermediate redshift galaxies in paper II 
that would be difficult to fit with another geometrical fig- 
ure. The efficacy of this approach is addressed in more 
detail below and in the Appendix (§A.3). 

3.3. Structural parameters 

Image structure is most commonly quantified via bulge- 
disk decomposition, yielding a bulge-to-total ratio, B/T. 
We refrain from this approach here, for reasons which we 
alluded to in the Introduction. For example, B/T param- 
eter may be poorly defined for asymmetric and compact 
galaxies. Irregularities in the surface brightness profiles, 
which can be caused by asymmetric structure, rings, or 
lenses, also cause problems for bulge-disk decompositions. 
While Kent showed that the concentration parameter cor- 
relates well with the bulge-to-total ratio, this holds only 
for objects with B/T < 0.63. At larger values of B/T, 
bulge-disk decomposition fails for several objects in Kent's 
sample, resulting in galaxies of type SO - Sa being given 
extremely high values of B /T. Bulge-disk decomposition 
also becomes unreliable when galaxy disks are fainter than 
the bulges. 

It is worth noting again that these problems mainly 
arise from older, one-dimensional methods of decompo- 
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sition. The newer two-dimensional decomposition tech- 
niques are clearly successful at reproducing the observed 
light profiles, wi th remarkably small residuals (Schad e 
et al. 1995, 1996; |Simard 1998] |Marleau fe Simard lggf ). 
Still, there are physical situations where bulge-disk de- 
composition techniques in general become problematic, 
namely where the astrophysical reality is more complex 
than simple bulge-disk models. Some galaxies have central 
condensations better des cribed by an expo nential profile 
rather than an r 1 / 4 -law (Wyse et al. 1997); many galax- 
ies have strong bi-symmetries, such as bars; virtually all 
galaxies have varying degrees of asymmetry due to star 
formation, dust, or large-scale gravitational perturbations 
and lopsidedness. All of these features represent details 
that decomposition into bulge and disk components do not 
address correctly. Simple disk and bulge decomposition is 
also inadequate for disk galaxies wher e the luminosity pro- 
file deviates from a pure exponential ( Freeman 1970| ), e.g. 
type I and type II disks. (Type I disk profiles have an 
added component which contributes to the light just out- 
side the bulge region; the surface brightness of a type II 
profile shows the opposite behavior (an inner truncation), 
and drops below the level of an exponential profile in the 
region near the center.) 

Given the astrophysical complexity of real galaxies, the 
physical interpretation of the derived model parameters 
of disk-bulge fits remains uncertain. Nonetheless, such 
profile-fitting methods should be useful for estimating non- 
parametric structural and photometric parameters (e.g. 
characteristic sizes, surface brightness, image concentra- 
tion, and ellipticity) in a way that uses the data in an 
optimal manner. In the current effort, however, we have 
taken a completely non-parametric approach of measuring 
sizes, surface-brightness, image concentration and asym- 
metry using multi-aperture photometry rather than deriv- 
ing a model-dependent B/T parameter. 

3.3.1. Half-light radii and surface-brightness 

We define first our working definition of a total mag- 
nitude since it represents the critical zeropoint for mea- 
surement of the half-light radius and surface-brightness. 
We use the dimensionless parameter rj to define the total 
aperture of the galaxies - a limiting radius which is not 
based on isophotes. 2 The concept of defining the size of 
a galaxy based on the rate of change in the enclosed light 
as a function of radius was first introduced by Petrosian 
(1976). In terms of intensity, rj can be defined as the ratio 
of the average surface bright ness within radius r to the 



local surface brightness at r ( Djorgovski fc Spinrad 1981 



of the light profile. These properties make it advantageous 
for faint object photometry. We defined the "total" aper- 
ture of the intermediate-redshift objects as twice the ra- 
dius r{rj = 0.2). The apparent total magnitudes arc then 
defined within this aperture. For ideal Gaussian or ex- 
ponential profiles, the magnitude mo. 2 within the radius 
2r(ij — 0.2) is approximately equal to the true total mag- 
nitude m t ou more than 99% of the light is included with 
the radius r(rj = 0.2). For an r 1 / 4 -law profile, there is a 
difference m .2 — m to t ~ 0.13 mag; this is due to the slow 
decline in luminosity at large radii that characterizes this 
profile. The radius r(rj — 0.2) was chosen based on vi- 
sual inspection of the curves of growth, derived from the 
aperture photometry, out to large radii. 

For reference, the theoretical value for the ratio of 
r(r) = 0.2) to half-light radius is 2.16, 1.95, and 1.82 
for three standard profiles: Exponential, Gaussian, and 
r 1 / 4 -law, respectively. The observed ratio is 2.3 ± 0.3 for 
B — V< 0.85 (with little trend with color), but rises slightly 
(2.6 ± 0.25) for the reddest galaxies with B-V> 0.85. A 
contributing cause to this rise is that for about half of the 
reddest objects, n/2 has been underestimated by ~ 20% 
because of their higher ellipticity. As we show in the Ap- 
pendix (§ A. 3. 2), the half-light radii of early-type galaxies 
with axis ratio a/b > 2 are systematically underestimated 
by up to 30%. This effect will also cause small changes to 
the measured image concentration (§ 3.3.2) of these galax- 
ies. 

A weak downward trend can be seen from blue towards 
red colors; this is what we expect since bluer objects tend 
to have exponential luminosity profiles, and redder objects 
are better described by r 1 / 4 -law profiles. However, this 
trend is broken by the reddest objects (B — V > 0.87), 
which have higher values of r(rj = 0.2)/rw2 than what is 

expected for an r 1//4 -law profile. 

Finally, the angular half-light radii r e were determined 
from the normalized curves of growth. Based on Mb and 
(corrected) R e we calculated the photometric-structural 
parameter SB e , the average blue surface brightness within 
the half-light radius, for all objects. For the nearby galaxy 
sample, the Tully catalog distances (as described in § 3.1) 
were used to determine R e (kpc). 

3.3.2. Image concentration 

We use the image concentration parameter C as defined 
by Kent (1985), which is based on the curve of growth. 
This parameter was shown to be closely correlated with 
Hubble type for "normal" galaxies: 



Sandage fc Pcrelmutcr 1990|) . Like Wirth et al. (1994), we 



C = 5log(r /n) 

In the above equation, r Q and fj are the outer and inner 
radii, enclosing some fraction of the total flux. In con- 
trast, the concentration parameter defined by Abraham 
et al. (1994) is not based on curve of growth radii, but on 
a flux ratio within two isophotal radii. 

However, in practice Kent also uses isophotes: He re- 
places the outer radius r a , which encloses 80% of the total 
light, by the radius of the 24th mag/arcsec 2 isophote. He 

2 Isophotal radii introduce redshift-dependent biases unless careful consideration and corrections are made for dimming due to the expansion 
(oc (1 + z)~ 3 in broad-band photon counts) and fe-corrections. While such redshift-dependent biases are not an issue for the samples studied 
in this paper, in future papers in the series this would be an issue were we not to avoid isophotes. 



follow Kron's (1995) suggestion to use the inverted form, 
r](r) = I(r)/(I(r)), which equals one at the center of the 
galaxy and approaches zero at large galactic radii. The 
radius r(ry = 0.5) corresponds roughly to the half-light 
radius r e . 

Since rj is defined as an intensity ratio, it is not affected 
by the surface brightness dimming effect that makes the 
use of isophotes problematic. Moreover, r\ is only depen- 
dent on the surface brightness within a given radius and 
not on any prior knowledge of total luminosity or the shape 
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has demonstrated that this radius encloses ~ 79% of the 
total light for all galaxy types in the restrictive confines 
of his sample ( Kent 1984 ) . Because of the surface bright- 
ness dimming effect that becomes important for non-local 
galaxies, we instead use a method that is independent of 
isophotes. The total aperture of the galaxy, which deter- 
mines the curve of growth, is defined based on the ^-radius 
as described in § 3.3.1. 

We have also explored the possibility of using 77-radii 
to define a concentration parameter. However, a concen- 
tration parameter based on the curve of growth was ulti- 
mately found to be the more robust measure: the curve of 
growth increases monotonically with galactic radius for all 
objects, while the 7y(r)-function will be non-monotonic for 
a "bumpy" light profile (like that of a well-resolved spiral 
galaxy). As a consequence, image concentration defined 
by the curve of growth rather than 77 exhibits less scat- 
ter when plotted against other correlated observables (e.g. 
color, surface-brightness) than an image concentration pa- 
rameter based on the 77-function. 

Anticipating our need to measure image concentration 
for small galaxies in paper II and future papers in this se- 
ries, we have studied the effects of spatial resolution and 
S/N on C. Here we focus primarily on resolution, as this 
was the dominant effect. The importance of resolution is 
demonstrated by the comparison of Schade et al. (1996) 
of decompositions of compact objects in ground-based and 
HST images: the cores of the blue nucleated galaxies are 
not resolved in ground-based imaging, and hence they are 
frequently misclassified as having much lower B/T-ratios 
than what is revealed by TJST-imaging. In paper II we 
analyze this sample of galaxies, and hence this illustration 
is of particular relevance. 

Resolution effects on image concentration were esti- 
mated by block-averaging the images of nearby galaxy 
sample over a range of values until the spatial sampling 
(as measured in pixels per half-light radius) was compara- 
ble to that of the compact galaxies at intermediate redshift 
observed with the WFPC2. The details of these simula- 
tions are presented in the Appendix (§A.l). In short, as 
the objects' half-light radii get smaller, the scatter in the 
measured concentration indices increases. While larger 
inner radii or a smaller outer radii decrease this scatter 
(due to improved resolution and S/N , respectively), such 
choices decrease the dynamic range of the concentration 
index. 

Based on these simulations, we chose a definition of 
C that is, to first order, sufficiently robust to allow a 
direct comparison of the image concentration of the lo- 
cal and the higher-redshift samples studied here and in 
paper II, and furthermore gives a large dynamic range: 
C = 5log(r(80%)/r(20%)). This concentration index is 
remarkably stable: The mean concentration does not de- 
viate from that measured in the original image by more 
than 0.2, or ~ 8% of the dynamic range in C, down to 
resolution of five pixels per half-light radius. 

Our definition is sufficiently close to that of Kent's 
(1985) so that it is meaningful to compare our values di- 
rectly to those he determined from photometric analysis 
of a sample of nearby galaxies. With this choice of radii, 
a theoretical r 1 / 4 -law profile has C — 5.2, an exponential 
profile has C — 2.7, and a Gaussian has C = 2.1. These 
values agree well with the results of Kent's analysis: he 



finds that elliptical galaxies have C ~ 5.2, and late-type 
spirals have C ~ 3.3. 

Lastly, since we use circular apertures, the measured im- 
age concentration may be affected by the ellipticity of the 
galaxy. Based on the comparison between our results for 
the Frei et al. sample and those of Kent's elliptical aper- 
ture photometry, we believe this to be a negligible effect 
in all cases but the earliest, must elliptic galaxies. Wirth, 
Koo and Kron (1994) found that for an r 1 / 4 law profile 
with axis ratio b/a — 0.2, the change in C is less than 
5%. The effect appears to be larger in our study. A more 
detailed description of this possible systematic is given in 
the Appendix. 

3.3.3. Image asymmetry 

The last image structure parameter is rotational asym- 
metry, A, as defined by Conselice, Bershady & Jangren 
(2000). This definition differs from earlier methods in that 
the asymmetry is determined within a constant ry-radius of 
77 = 0.2, a noise correction is applied, and an iterative pro- 
cedure which minimizes A is used to define the center of ro- 
tation. This algorithm was tested to be robust to changes 
in spatial resolution and signal-to-noise by Conselice et al. 
(1999) using simulations similar to those described here 
for the concentration parameter C; the systematics with 
resolution are below 10% of the original value for galaxies 
in paper I and II here. 

3.3.4. Morphological k— corrections 

To obviate the issue of 'morphological' k-corrections, im- 
age structural parameters should ideally be measured at 
the same rest-frame wavelength for all objects. Antici- 
pating our needs to derive the structural parameters for 
intermediate-redshift objects in paper II (and future pa- 
pers in this series), we have adopted the following protocol: 
(i) For the nearby galaxy sample we use the images in the 
Bj and g bands to derive the primary local image structure 
parameters. (The rest-frame wavelengths sampled by the 
R, r band images correspond to bands redshifted into the 
near-infrared for the intermediate-redshift galaxies.) (ii) 
We use the multi-band images of the Frei et al. sample to 
determine corrections to compensate for the wavelength 
dependence of asymmetry, concentration, and half-light 
radius - as described in the Appendix (§ A. 2). For exam- 
ple, HST Wide Field Planetary Camera-2 (WFPC-2) im- 
ages in the isi4 band of objects between 0.3 < z < 0.8 cor- 
respond to first order to the rest-wavelength range of the 
Bj and g bands. Nonetheless, the effective rest-wavelength 
for such intermediate-redshift galaxies is typically slightly 
redward of rest-frame Bj and g bands. The corrections 
in § A. 2 are suitable for such samples, as well as higher 
redshift samples imaged in redder bands. 

4. RESULTS 

4.1. Mean Properties, Distributions and Correlations 

The mean properties for our six parameters (Mb, B — V, 
R e , SB e , C, and A) are listed in Table 3, as a function of 
Hubble Type. While we would like to move away from us- 
ing 'Hubble Types,' they are so ingrained in the astronom- 
ical culture that they are a useful point of departure. For 
clarity in the following discussion, we group these types to- 
gether into "Early" (E-S0), "Intermediate" (Sa-Sb), and 
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"Late" (Sc-Irr). These names are potentially misleading, 
of course, and so we encourage the reader to treat them as 
labels which evoke, at best, a well-conceived galaxy type, 
but not necessarily an evolutionary state. Clearly further 
sub-division could be made, but our current purposes are 
illustrative, not definitive. 

A typical approach to exploring the correlations in (and 
dimensionality of) a multivariate parameter space is prin- 
cipal component analysis. While this is valuable, it is not 
particularly instructive for a first understanding of the dis- 
tribution of different types of objects in the parameter 
space. We are interested both in correlations between ob- 
servables and in trends as a function of the qualitative 
Hubble-type. These correlations and trends need not be 
one and the same. For example, two observables can be 
uncorrelated but still exhibit a distribution segregated by 
Hubble type. To develop such an understanding, we there- 
fore inspected the 15 possible 2-dimensional projections of 
our 6-dimensional parameter space. 

To distill this information further, we considered that 
there are in fact three types of physically-distinct param- 
eters: 

1. spectral index (color): this parameter is purely pho- 
tometric, by which we mean there is no information 
about the shape of the light profile. There is also no 
scale information, i.e. the amplitude and size of the 
light profile is also unimportant. In the balance of 
this paper we will use "color" and "spectral index" 
interchangeably. 

2. form (A, C): these parameters are purely structural, 
by which we mean that they do not depend - to first 
order - on the amplitude or the shape of the spectral 
energy distribution, nor on the physical scale of the 
light distribution; they reflect only the shape of the 
light profile. 3 

3. scale (R e , and SB e ): these parameters are physi- 
cally distinct. Luminosity is purely photometric (by 
our above definition). Size, which we also refer to as 
a structural parameter, is influenced by image shape, 
i.e., depending on the definition of size, two galaxies 
with different light profile shapes can have relatively 
different sizes (see §3.3.1, for example). Surface- 
brightness is a hybrid, photometric and structural 
parameter; it is a function of size and luminosity. 
While surface-brightness is a ratio of luminosity to 
surface area, it is still a measure of "scale" - in this 
case, the luminosity surface-density. 4 

This reduces the types of combinations (by parameter- 
type) to 6, i.e. between color, form, and scale. 

We find the strongest and physically most interesting 
correlations are between color, form, and the one scale pa- 
rameter, SB e (Figures 3-5). We focus on these for the 
remainder of the paper. Before turning to them, for com- 
pleteness we first summarize our observations of the other 
types of correlations: 

3 We consider image concentration to be a form, in contrast to Morgan who used it as a surrogate for spectral index. 

4 A fourth scale parameter which we do not consider here is line-width, or some measure of the amplitude of the internal dynamics. 

5 A few of the local galaxies have values of C that are lower than the theoretical concentration for an exponential disk (the errors in C are 
< 0.02 for all of them). The majority of these objects are late-type spiral galaxies with prominent, bright regions of star formation in the spiral 
arms. The star-forming regions cause the image profiles to become less centrally concentrated than a simple disk profile. 



Color-color correlations are strong and well known (e.g. 
Figure 2). Effectively they add higher-order information 
about spectral type. Here we consider only B — V as a 
simple spectral index which effectively represents the first- 
order information of spectral type. In general, one might 
adopt several spectral indices, e.g. U — V and V — K, or 
a single index based on multi-colors. 

Color-scale correlations also have been explored in de- 
tail elsewhere, e.g. color-luminosity relationships, known 
to exist for all galaxy types in both the optical and near- 
infrared (Huchra, 1977b; Mobasher et al. 1986; Bershady 
1995). The limited dynamic range of the Frei et al. sample 
in size and luminosity (they are mostly large and luminous 
systems) preclude useful results being drawn here in this 
regard. For example, the correlation of color with size in 
this sample is subtle and depends in detail on how size is 
defined, as noted above. Form-scale correlations including 
size and luminosity are also difficult to assess for this sam- 
ple for the same reasons of limited dynamic range in scale. 
However, scale versus scale is an interesting diagnostic be- 
cause, for example, size and luminosity allow one to probe 
the range of surface-brightness in the sample. We explore 
this in paper II. 

4.1.1. Spectral index versus form and scale 

Strong correlations exist in all three plots of color ver- 
sus form parameters C and A and scale parameter SB e 
(Figure 3). Early-type galaxies are redder, more concen- 
trated, high-surface-brightness, and more symmetric than 
Intermediate- and Late-type systems. The best correlation 
is between color and concentration in the sense that there 
is a smooth change in both quantities with Hubble Type. 
This is expected from a simple interpretation of the Hub- 
ble Sequence as a sequence parameterized by the relative 
dominance of a red, concentrated bulge (or spheroid) ver- 
sus a bluer, more diffuse disk. 5 In contrast, the distinction 
between Hubble Types in SB e and A is most pronounced 
between Early- types and the remainder; Intermediate- and 
Late-types galaxies are not well distinguished by either of 
these parameters. 

A more complete local sample will likely include a larger 
fraction of objects that do not follow these trends. For ex- 
ample, amorphous galaxies have surface brightnesses com- 
parable t o elliptical galaxies but a r e generally quite blu e 
in color ( |Gallagher fc Hunter 1987| ; pVIarlowe et al. 1997| ). 



Nonetheless, what is physically compelling about these 
color-form correlations is that each axis carries distinct 
information, respectively, on the integrated stellar popu- 
lation and its spatial distribution. 

4.1.2. Form versus form and scale 

There are clear trends present in the two plots of form 
versus SB e (scale) in Figure 4 as well as the plot of form 
parameters along in Figure 5. More centrally concen- 
trated galaxies have higher average surface-brightnesses 
and lower asymmetry; more symmetric objects have higher 
surface-brightness. In general, the concentrated, high 
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surface-brightness galaxies are Early-type, while the Late- 
type galaxies are less-concentrated, have lower surface- 
brightness, and are more asymmetric. While there is sub- 
stantial scatter in the form and scale parameters for Early 
and Late types, these two extreme groups still are well- 
separated in the above three plots. The Intermediate-type 
galaxies, however, are not well separated from these ex- 
tremes, and tend to overlap substantially with the Late- 
type galaxies, consistent with what is found in plots of 
color versus form and scale: Intermediate- and late-type 
galaxies have comparable degrees of asymmetry, and sim- 
ilar surface-brightness. 

One should be cautious in concluding the relative mer- 
its of form-scale and form-form and versus color-form and 
color-scale correlations based on the relative separation of 
Hubble Types. Using Hubble Types may be unfair if, for 
example, they were designed to correlate well with color 
but not necessarily with the quantitative form and scale 
parameters explored here. Since the form- form and form- 
scale correlations themselves are comparable, and nearly 
as strong as for color-form and color-scale, we are inclined 
to consider both as part of a general classification scheme. 
Certainly the form and scale parameters will each have 
different sensitivity to stellar evolution than color and so 
are advantageous to consider in isolation. 

4.1.3. Comparisons to previous work 

The correlation between image concentration and mean 
surface-brightness within the effective radius (Figure 4) 
has been explored by several groups in the context of 



galaxy classification ( Okamura et 
198l'Doi et al. 19931 



I. 1984 



Abraham et al. 1994 ) . We focus here 



Watanabe et al\ 



however, on Kent's (1985) r-band study since his definition 
of image concentration and effective surface-brightness are 
the most similar to our own. While similar, nonetheless 
the slope of the correlation is steeper for our sample, al- 
beit with much larger scatter, as illustrated in the top 
panel of Figure 6. As the middle and bottom panels re- 
veal, the cause of the steeper slope in our sample is due to 
a smaller dynamic range in image concentration. This is 
likely due to the fact that we use circular apertures when 
performing surface photometry, whereas Kent used ellip- 
tical apertures. We attempt to quantify the systematics 
due to differences in aperture shape in §A.3. While the 
dynamic range in image concentration is reduced using 
circular apertures for the Frei et al. sample, there does ap- 
pear to be a somewhat smaller scatter in C as a function 
of B - V. 

The nature of the large scatter in the top two panels of 
Figure 6 for the Frei et al. sample is also discussed further 
in §A.3. In short, we believe much of this scatter is due to 
uncertainties in the R- and r-band zeropoints of the Frei 
et al. sample. These uncertainties adversely affect only 
the surface-brightness values in Figure 6. Robust estima- 
tors of the scatter about a mean regression (i.e., iterative, 
sigma-clipping of outlying points) eliminate the outlying 
points, but still yield 50% larger scatter in i?-band SB e for 
the Frei et al. sample as a function of either image concen- 
tration or B — V. A plausible additional source contribut- 
ing to this larger scatter is that Kent's observed surface 
brightnesses were converted to face-on values, while ours 
were not "corrected" in this way. We conclude that if ac- 
curate and appropriate inclination corrections are possible 



to apply to all galaxies in a given study, this would be de- 
sirable. Since such corrections cannot be performed for the 
intcrmediate-redshift objects in paper II (and in general, 
if such corrections are not possible for a critical subset of 
the data), we believe it is best not make such corrections 
at any redshift. 

The asymmetry-concentration plane has also been ex- 
plored for galaxy classification purposes by, e.g., Abraham 
et al. (1994, 1996a) and Brinchmann et al. (1998). Our 
methods of measuring these parameters differ from theirs, 
and thus our quantitative results cannot be directly com- 
pared. However, a qualitative comparison to the A—C plot 
of Brinchmann et al. shows that both methods yield very 
similar results: the distribution of galaxies can be sub- 
divided into sectors where early-type, intermediate-type, 
and late-type objects dominate. Brinchmann et al. also 
use the local sample from the Frei catalog to define these 
bins, but note however that the points they plot represent 
a sample of intermediate-redshift galaxies. The A — C cor- 
relation in the Brinchmann et al. diagram is not as clear 
as that seen here for the local sample in Figure 5; the scat- 
ter in their diagram is comparable to the dynamic range 
of the parameters. This is probably due to the different 
properties of the samples, rather than to the differences 
in how we determine the parameters. For a more direct 
comparison, we plot B band asymmetry and concentration 
versus rest-frame B — V color for 70 galaxies from the Frei 
et al. sample (Figure 7), using both the A, C values from 
this study and those found by Brinchmann et al. It can be 
seen that the distributions are overall quite similar; how- 
ever, the separation in asymmetry of the different Hubble 
types is more apparent in this study, and the scatter in 
concentration is somewhat smaller. The conclusion here, 
then, is that our methodology offers typically modest, but 
sometimes significant improvements over previous work. 

4.2. Classification 

The above results point to how we can most effec- 
tively define a parameter demarcation to isolate, iden- 
tify, and classify normal galaxies. In the four-dimensional 
parameter space of (B — V, A, SB e , C), we define 
boundaries ("cuts") in the 6 two-dimensional projections 
between galaxies classified in the Hubble Sequence as 
Early/Intermediate and Intermediate/Late. These bound- 
aries, selected by eye on the basis of the distribution of 
Hubble types, are listed in Table 4 and illustrated in 
Figures 3-5. Segregation by higher-dimensional hyper- 
surfaces are likely to be more effective (galaxies appear 
to be distributed on a 'fundamental' hyper-surface - the 
subject of a future paper), but the projected boundaries 
here are meant as illustrative, and practical for application 
when all of the parameters are not available. We stress 
that these boundaries are not definitive in some deeper 
physical sense. For example, in terms of formal Hubble 
Types cuts involving color are clearly "best;" however, as 
noted above, this may not be physically significant. 

It would be uninteresting if all of the cuts provided the 
same classification. Moreover, one expects there will be 
discrepancies for objects near boundaries. We find that 
49% of the sample matches in all cuts, while 64%, 87%, 
and 99% of the sample matches in at least 5, 4, or 3 cuts, 
respectively. (Hereafter, we refer to cases where 5 out of 
6 cuts match as "5/6," etc.) This degree of consistency 
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seems reasonable so we have not tried to fine-tune the 
boundaries (such fine-tuning would not be sensible any- 
way since the details of the classification self-consistency 
are likely to be sample-dependent): The preponderance of 
objects are classifiable by a simple majority of the classi- 
fications based on the 6 cuts; 13% of the objects have a 
more ambiguous classification. 

Of interest are the discrepancies within and between 
cuts in different combinations of color, form, and scale. 
We found that it is useful to group the six cuts into two 
groups of three. The first consists of the cuts in Fig- 
ure 3 between color, form and scale, which we refer to 
as color- form/scale. The second consists of the cuts in 
Figures 4 and 5 between form and scale, which we refer 
to as scale/form- form. For example, 64% of the variance 
in the 5/6 cases comes from cuts in C-SB e , whereas cuts 
in C (B — V) and A-C are always consistent with the 
majority classification. More generally, scale/form- form 
cuts are internally mis-matched 40% of the time, while 
color-form/scale cuts are internally mis- matched only 21% 
of the time (and two-thirds of these color- form/scale mis- 
matches are also present in scale/form- form mis- matches). 
In other words, the color- form/scale cuts tend to be more 
consistent; much of the variance in the scale/form- form 
cuts again comes from C-SB e . 

Only two galaxies pose a substantial problem for classi- 
fication: NGC 4013 and NGC 4216. They are classified by 
various cuts to be in all categories (Early, Intermediate, 
and Late), and have no majority classification. However, 
both arc highly inclined (4013 is edge on), which appears 
to give them unusual observed properties. Indeed, they 
are extreme outliers in several of the projections in Fig- 
ures 3 and 4 (see also A. 3. 2 and figures therein). Hence 
such problem cases are likely to be easy to identify. Three 
other sources classified in all three categories (NGC 4414, 
4651, and 5033) are not a problem: They have 4/6 con- 
sistent classifications. Two of these (NGC 4651 and 5033) 
have Seyfert nuclei, and are outliers only in plots with im- 
age concentration; they are highly concentrated for their 
color. NGC 4414 is not an outlier in any of the plots. 

Finally, it is interesting to note that 23% of the 
sample has inconsistent majority classifications in color- 
form/scale versus scale/form- form cuts. This is true for 
100% of the 3/6 cases, and 55% of the 4/6 cases. However, 
we believe this is for different reasons. In the latter cases 
(only) we find that the galaxies are predominantly at high 
inclination (^50% excess in the top half and top quartile 
of the sample distribution in inclination). Moreover, the 
color-form/scale classifications in these cases are all ear- 
lier than the majority scale/form- form classifications. We 
surmise this is due to the effects of reddening on B — V . 6 
While the color- form/scale classifications tend to be earlier 
for the 3/6 cases, because there is no apparent inclination 
dependence, these differences are due likely to other phys- 
ical effects. Two possibilities include low star-formation 
rates or high mctallicity for galaxies of their form. Both 
of these conjectures are testable via spectroscopic obser- 
vation. 

We suggest then, as a practical, simple prescription, 



that the majority classification for all 6 cuts be taken as 
the classifier, except in the situation where the galaxy in 
question is highly inclined. In the latter case, the ma- 
jority classification of the scale/form- form cuts should be 
adopted. When galaxies have only 3/3 consistent clas- 
sifications, (13% of the Frei et al. sample), the adopted 
classifier should be intermediate between the two most 
common classifications. It also may be of interest to note 
if the color- form/scale and scale/form- form majority clas- 
sifications differ. However, further elaboration based on 
these two-dimensional projections of a higher-dimensional 
distribution is not likely to be warranted. 

4.2.1. Discussion 

We note that there are no distance-dependent scale pa- 
rameters in our classification. By this we mean specifically 
that the classification parameters do not depend on knowl- 
edge of the distance modulus. Hence this classification 
is both quantitative and independent of the cosmological 
distance-scale and its change with cosmological epoch (i.e. 
no a priori knowledge is needed about Ho or qo). The ef- 
fects of the expansion do change the observed classification 
parameters. However, with knowledge of galaxy redshifts 
and judicious choice of "rcdshifted" photometric bands, 
surface-brightness dimming can be corrected and band- 
shifting cither eliminated or corrected via the protocol de- 
scribed in the Appendix. Galaxy evolution, of course, will 
also modify the values of the parameters, but this is pre- 
cisely the utility of the classification systems as applied to 
such a study: In what way do the parameters and their 
correlations evolve? How do the scale parameters change 
for a fixed range in classification parameters? These are 
issues which we intend to explore in subsequent papers in 
this series. 

We also comment on the efficacy of using the four- 
dimensional parameter space of color, concentration, 
surface-brightness, and asymmetry for the classification of 
distant galaxies. As noted earlier, Abraham et al. (1996a) 
and Brinchmann et al. (1998) have explored the use of 
the asymmetry-concentration plane as a tool for distant 
galaxy classification. The use of the additional parameters 
of color and surface-brightness are clearly advantageous; 
they offer substantially more information, particularly as 
a diagnostic of the stellar population age and surface- 
density. The reasoning behind using A and C alone has 
been that to first order, they can be estimated without red- 
shift information. Yet the wavelength dependence of both 
parameters (i.e., what is referred to as 'morphological k- 
corrections') can lead to measurement systematics. These 
systematics, if not corrected, in turn result in objects over 
a range in redshift to be systematically misclassified. For 
example, Brinchmann estimates that at z = 0.9, 25% of 
spiral galaxies are mzs-classified as peculiar objects in the 
A — C plane. This fraction is expected to increase at larger 
redshifts. Hence, for high-z studies of galaxy morphology, 
redshift information is crucial even when using asymme- 
try and concentration. Therefore, since redshift informa- 
tion is crucial no matter what, there is no reason not to 
use the four-dimensional classification we have outlined in 



inclination will also cause changes in other measured parameters. Changes in C, however, appear to be small (see §A.3.2). Surface- 
brightness will tend to increase at modest inclinations, and then decrease at high inclinations if a prominent dust lane obscures the bulge. 
Likewise, A may increase due to a dust lane until the galaxy is directly edge-on. As a consequence of these changes and the distributions and 
cuts, C-SB e tends to mimic the color-form cuts in the high-inclination cases, while A-SB e and A-C do not. 
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future studies. The recent refinements and calibration of 
the technique of estimating redshifts photometrically make 
this all the more tractable. 

Finally, we note that while the classification we have 
proposed here is practical and useful, there are five areas 
where we anticipate it can be improved or elaborated: (i) 
As we have mentioned before, the spectral-index parame- 
ter could have much greater leverage at distinguishing be- 
tween different stellar populations by adding pass-bands 
that expand the wavelength baseline (e.g. the U and K 
bands in the near-UV and near-IR, respectively), or by 
increasing the spectral resolution (e.g. line-strengths and 
ratios). A further step of elaboration would be to explore 
spatially resolved spectral indices (gradients) and deter- 
mine their correlation with form parameters, (ii) Internal 
kinematics should be considered. Ideally, the kinematic 
information would include estimates of both the random 
and ordered motion (rotation) so that the dynamical tem- 
perature could be assessed, in addition to the overall scale. 
Kinematics are relatively expensive to obtain (compared 
to images), but with modern spectrographs on large tele- 
scopes, the absolute cost is minimal at least for nearby 
galaxies, (iii) Higher-dimensional correlations are worthy 
of exploration to determine, for example, whether "fun- 
damental" hyper-planes can adequately describe the en- 
tirety of the galaxy distribution, (iv) It is worth consider- 
ing whether there are additional form parameters of value 
for classification that have not been included here, (v) 
The classification scheme needs to be tested against much 
larger, and more volume-representative samples of galax- 
ies. 

5. SUMMARY AND CONCLUSIONS 

We have presented results from a study of the photomet- 
ric and image-structural characteristics and correlations of 
a sample of local, bright galaxies (Frei et al. 1996). We find 
it illuminating to distinguish between parameters which 
characterize spectral-index (color), form (image concen- 
tration and asymmetry), and scale (size, luminosity, and 
surface-brightness). In this context, we arrive at the fol- 
lowing main results and conclusions. 

• We find that a combination of spectral-index, form 
and scale parameters has the greatest discrimina- 
tory power in separating normal Hubble-types. The 
strongest correlation is found between color and im- 
age concentration. However, there are equally strong 
correlations between form parameters (e.g. A and 
C), but here the Hubble-types are not as well dis- 
tinguished. As an indicator of classification utility, 
we suggest that the strength of the correlation be- 
tween parameters is likely more important than the 
separation of Hubble Types within the correlation. 



• It is possible to define a quantitative classifica- 
tion system for normal galaxies based on a four- 
parameter sub-set of spectral-index, form and scale: 
rest-frame B — V color, image concentration, asym- 
metry, and average surface brightness within the 
half-light radius. We propose a specific classifica- 
tion that distinguishes between "normal" galaxies 
as Early, Intermediate, and Late based on cuts in 
these four parameters. The classification is success- 
ful for 99% of the Frei et al. sample. Nonetheless, 
we designate this as "preliminary" until larger, more 
comprehensive samples of galaxies are needed than 
analyzed in the present study. 

• Distance-dependent scale parameters are not part of 
this preliminary classification. 

• These classification parameters can be measured re- 
liably over a broad range in S/N and image resolu- 
tion, and hence should be applicable to reliably dis- 
tinguishing between a wide variety of galaxies over 
a large range in redshift. 

• Redshift information is needed to estimate reliably 
both the photometric properties (rest-frame color 
and surface brightness) as well as the structural pa- 
rameters asymmetry and concentration at a fixed 
(B-band) rest-frame wavelength. In terms of redshift 
independence, asymmetry and concentration alone 
thus offer no advantages over the additional param- 
eters classifiers proposed here. Indeed, incorporating 
the full suite of parameters defined here is advanta- 
geous for the purposes of classification. 
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APPENDIX 

A. CORRECTIONS FOR MEASUREMENT SYSTEMATICS 

Here we establish the measurement systematics due to changes in image resolution for half-light radius and image 
concentration, and for band-shifting effects on half-light radius and image concentration, and asymmetry. Systematic 
effects of images resolution and noise on asymmetry are quantified in Conselice et al. (1999). 



A.l. Resolution dependence of observed size and image concentration 
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To maximize the dynamic range of the measured concentration index, C, the inner radius should be small, and the 
outer radius large relative to the half-light radius. In this way, one samples the light profile gradients in both the central 
and outer regions of a galaxy where the bulge and disk contribute quite differently. This strategy maximizes the leverage 
for discriminating between different profiles, e.g. exponential and r 1 / 4 -law. In the presence of noise and limited spatial 
resolution, however, the choice of radii determines the robustness of the concentration index: As noted by Kent (1985), 
the inner radius should be large enough to be relatively insensitive to seeing effects, and the outer radius should not be so 
large that it is affected by uncertainties in the sky background and S/N. In the current study, the sources are resolved, 
and the images are at moderately high signal-to-noise: within the half-light radius, the sample of local galaxies have 
600 <; S/N <; 3000. The intermediate redshift galaxies in paper II have S/N in the range 40 to 90, with a mean of ~ 55. 
This is sufficiently high that we focus our attention here on the effects of spatial sampling and resolution. 

Even in the absence of significant image aberrations, an additional limiting factor is the number of resolution elements 
sampling the inner radius. This is likely to become a limiting factor when the half-light radii is only sampled by a few 
pixels. To understand this potential systematic, we have calculated six concentration indices C = 5log(r a /ri) for several 
different choices of inner and outer radii. We use enclosing 20 and 30% of the light, and r D enclosing 50, 70, an 80% of 
the light. The radii were measured for nearby galaxies that were block-averaged by factors 2, 4, and 6 to simulate coarser 
spatial sampling, as shown in Figure 8. The six different concentration indices are plotted as a function of sampling 
in Figure 9. These simulations span sufficient dynamic range in size to cover most galaxies observed, for example, in 
the Hubble Deep Field. With factors of 4 and 6, we measure radii with pixel sampling similar to that observed in the 
HST/WFPC-2 images of the intermediate-redshift objects of paper II. Typically, these galaxies have half-light radii of 
0.3 - 0.7". For the Planetary Camera, the scale is 0.046"/pixel, and for the Wide Field, 0.10"/pixel; hence the half-light 
radii are of order 3 to 15 pixels. 

The half-light radius r(50%) is remarkably stable, even with poor sampling. Unfortunately, the dynamic range given 
by concentration indices with r — r(50%) is too small to be useful. As expected, the 30% radius was more stable than 
the 20% radius to decreased spatial resolution. However, the concentration indices using r(30%) were less sensitive to 
the differences between galaxy types, and gave a smaller dynamic range than indices using r(20%). The inner radius 
dominated the effect on the amplitude of the systematics; changing the outer radius from 70% to 80% decreased the 
scatter only marginally. With a block-averaging factor of 6, where the half-light radii are typically only ~ 5 pixels, the 
scatter becomes large for all choices of concentration indices. 

Based on these simulations, we decided to use the radii enclosing 80% and 20% of the total light (as did Kent, 1985), even 
though r D = r(70%) gives concentration indices with slightly smaller scatter at poor resolution. For objects with half-light 
radii of only 7 pixels, the mean differences in concentration (relative to the original image) are ACso:20 = — O.lOig (so an( ^ 
AC70:2o = — O.IO^Q so- Even at a resolution of only five pixels per half-light radius, the concentration index only deviates 
by 0.2 relative to the original image; this is ~ 8% of the dynamic range in Cso : 20- Thus we consider this parameter to be 
robust enough to useful in the comparison of local and intermediate-redshift samples. 

A. 2. Systematics with wavelength 

Observations at different wavelengths sample preferentially different stellar populations in a galaxy. Since these popula- 
tions are not always spatially homogeneous, the image-structural characteristics (concentration, asymmetry, and half-light 
radius) will will have some wavelength dependence [e.g., see de Jong's (1995) study of disk scale- lengths] . Hence, when 
comparing one of these parameters for different galaxies, the parameter ideally should be measured at the same rest-frame 
wavelength for all objects. This is not possible in general for studies over a wide range in rcdshifts employing a finite 
number of observed bands. To determine the amplitude of the wavelength-dependence for the measured structural param- 
eters, we therefore compare the Bj and R structural parameters for 72 of the Frei et al. galaxies. The differences between 
the red and blue structural parameters versus the rest-frame color, B — V, arc shown in Figure 10. For comparison, all 
intermediate-redshift objects in paper II, except two, fall in the bluest bin (B — V < 0.62). 

The plot of AC = Cb — Cr shows that in most cases, the values are slightly negative, i.e. the majority of objects are 
more highly concentrated in the red band than in the blue, as expected because of the redness of the central bulge. Only 
the bluest galaxies have comparable image concentration in both bands. There is a weak trend towards more negative 
values for the redder (early-type) objects, which also show a larger scatter than the bluer objects. 

In the plot of AA = Ab — Ar it is clear that most galaxies have at positive values, i.e. their image structures are more 
asymmetric in the blue band than in the red, as shown by Conselice (1997). The difference in asymmetry is only seen 
for late- and intermediate-type objects; red objects are generally very symmetric in both bands, and have Ab — Ar ~ 0. 
This trend was also noted by Brinchmann et al. (1998). 

The plot of half-light radii (AR e = R £ _b — Re,n) shows that most values are slightly positive, with a larger scatter for 
redder objects. No other trend with color is seen. The fact that the galaxies have slightly larger half-light radii in the 
blue band is consistent with their image concentration being higher in the red band, as a bulge profile generally has a 
much smaller scale length than an exponential profile. 

In summary, the average differences (±lcr) between parameters for galaxies with B — V < 0.62, determined from the 
Bj and R bands, are: AC = -0.15 ± 0.30, AA = 0.013 ± 0.044, and AR e = 0.21 ± 0.80 kpc. 

A. 2.1. Corrections for wavelength systematics 

Based on the mean values above, we correct the measured structural parameters for galaxies at non-zero redshift 
to the rest-frame B band values as follows, where for clarity we use the intermediate redshift galaxies in paper II as 
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an example. The structural parameters of these intermediate-redshift objects generally were measured at rest-frame 
wavelengths between Bj and R, i.e., in the observed WFPC2 J 81 4 band for z ~ 0.6. For "normal" galaxies, this would 
cause us to overestimate C, and underestimate A and R e . Hence we use the differences listed above and the redshift of the 
objects to linearly interpolate the correction to the measured values. Specifically, for a given parameter and color bin, we 
use the mean difference between values measured in the Bj and R bands, and the position of the rest-frame wavelength 
relative to the Bj band, to make corrections to the measured parameters. (Note that the correction made to R e also 
affects the value of SB e in general.) For some objects, the combination of observed band-pass and redshift corresponds 
to rest- frame wavelengths slightly blueward of Bj. When computing the corrections for these objects, we assumed that 
the wavelength trends continue outside the Bj — R wavelength range. Overall, these corrections are small for objects in 
paper II, while for higher-redshift objects we expect band-shifting effects to become increasingly important. 

We add a final, cautionary note that it is not certain the corrections for intermediate-redshift objects should be made 
based on the correlations we see for the nearby sample. When comparing the observations in the bluer bands (-B450 or 
Vfioe) to those in, e.g. the Ism used in paper II, we find that most objects are more concentrated in the blue band, and 
slightly larger in the red band - this is the opposite of what we see for the Frei et al. sample.^] For asymmetry, the trend is 
the same for both samples (higher A in bluer bands). The trends are not directly comparable, however, to what we see in 
the local sample, as the observations in the bluer bands correspond to rest-frame wavelengths in the UV region for most 
intermediate-redshift galaxies. For this reason, and since the small sample of intermediate- z objects poorly defines the 
variation in image structure with wavelength, we adopt the more well-determined trends seen for the Frei et al. sample 
to calculate the band-shifting corrections. These corrections based on local galaxy trends tend to make the intermediate- 
redshift objects somewhat less "extreme" ; their half-light radii become larger, their surface brightnesses fainter, and their 
image concentrations lower. If instead we had based our corrections on the trends seen within the intermediate- z sample 
of paper II, then this sample would be even more extreme relative to the local galaxy sample. The corrections would then 
tend to shift the positions of the intermediate- z objects even farther from the nearby galaxies in diagrams that include 
any of the parameters R e , SB e , and C . 

A. 3. Systematics with aperture shape 
A. 3.1. Comparison to elliptical aperture photometry 

Circular-aperture surface-photometry will yield systematic differences in the measured structural parameters when 
compared to those derived from elliptical-aperture surface-photometry To assess this, we compared our results for the 
Frei et al. catalog in R and r bands to the results of Kent (1985) for a sample of local, Hubble-type galaxies (Figure 
6). Kent used elliptical apertures tailored to fit the axis ratio and position angle of each isophote in galaxy images to 
determine r band image concentration and average surface brightness within the half-light radius. 

As we detail in the figure caption, we have attempted to transform all of the surface-brightness values to the Cousins 
R band (R c ). For each of the relations in Figure 6 we have characterized the slopes and scatter about a mean regression 
using a simple linear, least-squares algorithm with an iterative, sigma-clipping routine to remove outlying points. Given 
the nature of the data, such an algorithm is not statistically correct (see, e.g., Akritas and Bershady, 1996). However, 
given the potential photometric uncertainties (discussed below) and the need for robust estimation, it is not possible 
to formally implement more appropriate algorithms. Nonetheless, the relative characterization of the slopes and scatter 
between Frei et al. and Kent samples is useful. 

As discussed in §4.1.3, the slope of the correlation between average surface-brightness and image concentration is steeper 
for our study than for Kent's because of a decreased range in image concentration in our study. The effect (bottom panel 
of Figure 6) is such that the bluest galaxies have comparable image concentration values in both studies while the image 
concentration of the reddest galaxies differ by as much as 1 unit in the mean (Kent's values are larger). We interpret 
this as likely to be the effect of different aperture shapes. The results of our study of systematics with axis ratio (below) 
support this conclusion. Surprisingly, there is no indication that elliptical apertures give significantly different results 
than circular apertures for intermediate- and late-type (disk dominated) galaxies. 

The larger scatter in the Frei et al. (1996) sample in the top two panels of Figure 6 might lead one to conclude that the 
elliptical apertures provide a superior measurement of effective surface-brightness. However, much of the scatter is due to 
the subset of the Frei et al. sample observed at Palomar Observatory. We believe that zeropoint problems are the cause 
of much of this scatter, consistent with discussion in Frei et al. concerning the difficulty of photometric calibration. The 
bulk of the objects observed at Lowell Observatory are consistent with independent i? c -band photometry from Buta and 
Williams (1996), although there are some points that are very discrepant. In general, the overlap is excellent in SB e and 
B — V between the Kent sample, the Lowell subset of the Frei et al. sample, and the subset of the Frei et al. sample with 
Buta and Williams' photometry. 

A.3.2. Systematics with axis ratio 

A second approach to determine the systematic effects of aperture shape on measured structural parameters was also 
used: we quantify the degree to which "normal" galaxies with the same intrinsic morphology but with different axial ratios 
a/b will have different C when measured with circular apertures. The galaxies in the Frei et al. catalog were divided into 

7 Indeed, Huchra noted that the Markarian galaxies get bluer toward their centers, reminiscent of the blue "bulges" seen in the blue nucleated 
galaxies of paper II, yet in contrast to the color gradients found for "normal" galaxies. This type of color-aperture relation was also noted by 
de Vaucouleurs (1960, 1961) for the latest Hubble-type galaxies (Sm, Im). 
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early-, intermediate- and late- type objects (using the same bins as elsewhere in the paper), and we plot concentration and 
half-light radius versus the logarithm of the axis ratio (taken from the RC3 catalog) . 

In the image concentration plots (Figure 11), a weak trend can be seen for the late- type galaxies (top panel), with 
slightly higher values of C for the more inclined objects. This effect, if caused by the shape of the apertures, will lead us 
to overestimate the image concentration by at most ~ 0.1 (3%) for the nearly edge-on galaxies. We do not expect this 
to be a problem for our analysis. The two labeled objects have unusually high values of C for their morphological type. 
One of them, NGC 5033, is known to be a Seyfert 1 galaxy; the other, NGC 4651, is a suspected "dwarf-Seyfert" galaxy 
(Ho et al. , 1997). For intermediate- type objects (middle panel), no trend is observed. The lowest C- value, which belongs 
to NGC 4013, could be caused by the prominent dust lane in this object: the central light distribution is divided into 
two parts, making it difficult to determine the position of the center. Effects like these will likely be more problematic 
for objects with high values of a/b. The highest C-value in this plot is that of NGC 4216, which also is highly inclined 
and has spiral arm dust lanes superimposed on the bulge. In the bottom panel, a trend is observed for the early-type 
galaxies: the concentration is lower for objects with higher a/b— ratio. This effect will cause us to underestimate the 
image concentration of these objects by ~ 0.5, or 10-15%. This result agrees well with what was seen in the comparison 
of the Frei et al. sample to Kent's image concentration measurements, as described above. This leads us to conclude that 
our circular aperture photometry will underestimate the image concentration somewhat for elliptical/SO galaxies. Again, 
there is no indication that the aperture shapes lead to different results for intermediate- and late-type galaxies. 

In the plots of half-light radius R e versus a/b (Figure 12), no trends are seen for the intermediate- and late-type objects. 
For the early-type objects, however, the measured half-light radii become progressively smaller for increasing values of 
a/b. The trend is weak; it will cause us to underestimate the half-light radii by at most 30 % for objects with a/b ~ 4. 
If this effect is real, the derived surface brightness will be too bright by < 0.7 magnitudes for the most highly elliptical 
early-type galaxies. 
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TABLE 1 

Nearby Galaxy Data Sets 



Survey 


N 


Distance 


Telescope 


Filters 


Frei et al. (1995) 


72 


< 41 Mpc 


Lowell 1.1m 


Bj, R 




29 


< 19 Mpc 


Palomar 1.3m 


gri 


Kent (1984) 


53 


< 65 Mpc 


Whipple 0.6m 


r 


Huchra (1977) 


196 


< 286 Mpc 


Mount Wilson 2.5m, 


UBV 








Palomar 1.5m 
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TABLE 2 

Local Galaxy Properties 
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,7±1. 


.2 


3 


,58±0 


,04 





,14±0 


.07 


4866 


-1 


-20 


17±0. 


.10 


0, 


,90±0. 


14 


20 


,35±0. 


.10 


2. 


,46±0. 


.03 


62. 


,7±0. 


.6 


3 


,91±0. 


,03 





,04±0 


.03 


5005 


4 


-21 


.86±0. 


.10 


0, 


,80±0. 


11 


20 


,57±0. 


.11 


5. 


,92±0 


.15 


96 


,7±0. 


.4 


3 


,91±0. 


.01 





.20±0 


.01 


oUoo 


r 



-21 


.54±0. 


,1U 


0, 


.55±0. 


15 


oo 

AA 


.U4±U. 


1 1 
,11 


1U 


,Uo±U. 




166 


,2±1. 


.4 


A 

4 


i n 

,oy±u. 


no 
.Uz 


u 


1 o i n 

,lo±0 


.Uo 


5055 


4 


-20 


.78±0. 


,10 


0, 


,72±0. 


14 


21 


,14±0. 


.10 


4 


.69±0 


.05 


216. 


.5±0. 


.6 


3. 


,45±0 


.00 





19±0 


.02 


5204 


9 


-17 


,58±0. 


,14 


0. 


,41±0. 


.18 


21 


,79±0. 


.15 


1. 


,45±0 


.04 


90. 


,5±1. 


.0 


2 


,88±0 


.07 





,25±0 


.06 


5248 


4 


-21 


,82±0. 


,13 


0. 


.65±0. 


.18 


21 


,34±0. 


.13 


8. 


,32±0 


.08 


103 


,6±0 


.3 


3 


,19±0. 


.02 





.20±0 


.02 


5322 


-5 


-22 


,34±0. 


,16 





,89±0. 


.18 


20 


,99±0. 


.18 


8. 


,97±0. 


.36 


105. 


.3±3. 


.2 


4, 


,87±0. 


,02 





,02±0 


.01 


5334 


5 


-21 


.01±0. 


.22 








23 


,41±0. 


.22 


14 


,84±0 


.23 


126 


,3±1. 


.4 


3 


,36±0 


.04 





,09±0 


.12 


5364 


4 


-21 


.84±0. 


.22 





,61±0. 


.18 


22 


.29±0 


.29 


12. 


,97±1 


.09 


161 


,5±0 


.7 


2 


.97±0 


.01 





15±0 


.02 


5371 


4 


-22 


.45±0. 


.14 


0. 


,68±0. 


.20 


21 


,64±0, 


.14 


12. 


,73±0. 


.06 


105 


.7±0. 


.2 


2. 


,61±0. 


,02 





,14±0 


.02 


5377 


1 


-21 


14±0. 


,18 


0. 


.89±0. 


.18 


21 


,38±0. 


.21 


6 


,18±0. 


.34 


76. 


,2±0 


.9 


4, 


,32±0 


,03 





.07±0 


.02 


5585 


7 


-18 


,79±0. 


,14 


0, 


,46±0. 


.20 


21 


.97±0. 


.15 


2. 


,74±0. 


.05 


135 


.8±0. 


.5 


2 


,80±0 


.02 





,16±0 


.03 
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TABLE 2-Continued 



NGC 


Type 


M B 


B-V 


SB e 


Re 


r v =o.2 


C 


A 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


5669 


6 


-20.85±0.21 




21.61±0.21 


6.00±0.06 


84.8±0.5 


2.74±0.04 


0.19±0.06 


5701 





-21.30±0.14 


0.84±0.20 


20.73±0.14 


4.91±0.04 


52.1±0.3 


4.28±0.01 


0.05±0.02 


5792 


3 


-21.35±0.13 


0.78±0.18 


22.19±0.13 


9.85±0.11 


90.1±0.7 


3.32±0.02 


0.36±0.03 


5813 


-5 


-21.74±0.15 


0.94±0.18 


21.02±0.17 


6.89±0.24 


83.4±0.9 


4.07±0.03 


0.03±0.02 


5850 


3 


-21.68±0.13 


0.74±0.18 


21.54±0.13 


8.51±0.09 


109.2±1.5 


3.99±0.04 


0.08±0.03 


5985 


3 


-22.09±0.13 


0.74±0.18 


22.20±0.13 


13.99±0.06 


107.1±0.5 


2.71±0.01 


0.11±0.02 


6015 


6 


-20.45±0.10 


0.56±0.14 


21.76±0.11 


5.34±0.07 


91.4±0.4 


2.70±0.01 


0.19±0.02 


6118 


6 


-21.10±0.13 


0.64±0.18 


22.38±0.13 


9.62±0.10 


112.8±0.7 


2.78±0.03 


0.14±0.04 


6384 


4 


-22.14±0.14 


0.61±0.20 


21.42±0.14 


9.95±0.06 


138.5±0.5 


3.27±0.02 


0.13±0.02 



Note. — Col. 1: the galaxy identification. The properties listed in cols. 5-9 are determined from the g and 
Bj band images. Col. 2: the morphological type of the galaxy. Col. 3: the absolute B band magnitude. Col. 
4: the rest-frame B — V color. Col. 5: the average surface brightness within the half-light radius. Col. 6: the 
half-light radius (in kpc). Col. 7: the radius corresponding to rj = 0.2 (in arcseconds). Col. 8: the image 
concentration. Col. 9: the image asymmetry. Quantities in columns 7-9 are for the observed B or g bands. 



TABLE 3 

Average Galaxy Properties 



Type 


Sample 


M B 


B - V 


Re 


SB e 


c 


A 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


E-SO 


Frei 


-2i-4±S:2 


o.9l8:g 


4.6±l:i 


2o.3i8;l 


4-218:3 


o.o3i8;8^ 


Sa-Sb 


Frei 


-2i 5+°- 6 


o.7l2:i 


s.3±i:2 


91 0+O.8 
^ 1 -°-0.2 


o r + 0.4 
°.°-0.5 


1S+ 003 


Sc-Irr 


Frei 


-2o.6±J:2 


o.6l2:i 


e.4±i:i 


2i.7±2;5 


o q+0.3 
°" U -0.3 


17+ - 07 

u - 1 '-0.04 


E-SO 


Kent 


-20.81^ 


o.9l2;2 






r q+0.3 




Sa-Sb 


Kent 




o.8l2:l 






4 2+0-5 




Sc-Irr 


Kent 


-20.91°;? 


o.6l2:l 






o 4 +0.4 




E-SO 


Huchra 


-20.5±g:I 


u -'-o.o 










Sa-Sb 


Huchra 


-21.01^ 


o-7l2:J 










Sc-Irr 


Huchra 


-21.0±S;i 


o.5l2:l 











Note. — Col. 1: the Hubble Type of the objects. Col. 2: a code representing 
the sample from which each source was extracted: Frei et al. (1995), Kent (1984), 
or Huchra (1977). Cols. 3-8: the same as as columns 3-8 in table 2. Structural 
parameters are corrected to correspond to the rest-frame _B-band except for Kent's 
sample where r band image concentration is listed. For all columns, the middle 
50% of the sample is contained within the listed bounds. 
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TABLE 4 

Classification Boundaries 



Parameters 


Early /Intermediate 


Intermediate/Late 


Figure 


(1) 


(2) 


(3) 


(4) 


SB e , B-V 


B-V = 0.12 SB e - 1.66 


B-V = 0.64 


4 


C,B-V 


B-V = -0.23 C+ 1.78 


B-V = -0.04 C + 0.76 


4 


A, B-V 


B-V = 1.05logA + 2.21 


B — V = 0.32 log A + 0.88 


4 


C, SB e 


SB e = 7.04 C - 7.02 


SB e = 3.52 C+ 10.73 


5 


A, SB e 


SB e = -10.4/o#,4 + 7.29 


SB e = -2.46/o 5 .4 + 19.36 


5 


A, C 


C = 21. 5 log A + 31.2 


C = 2.44 log A + 5.49 


6 



Note. — Col. 1: Parameter pair. Col. 2: the boundary between early and 
intermediate-type galaxies. Col. 3: the boundary between intermediate and late- 
type galaxies. Col. 4: the corresponding figure in this paper. 
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Figure captions 

Fig. 1. — Rest-frame B — V versus Mb for the nearby galaxy sample of Frei et al. (E-SO, Sa-Sb, and Sc-Irr) and Huchra's 
(1977a) sample of normal Markarian galaxies (pluses). The dotted outline indicates the approximate locus of dE/dSph 
galaxies. The intermediate redshift samples from paper II are also plotted for comparison: blue nucleated galaxies (BNGs) 
compact, narrow emssion-line galaxies (CNELGs), and small, blue galaxies (SBGs). (The two SBGs and the BNG that 
we ultimately determine not to be "Luminous Blue Compact Galaxies" in paper II are shown as hatched symbols.) Only 
a few Markarian galaxies and late-type galaxies from the Frei et al. catalog share the extreme color-magnitude properties 
of the intermediate-redshift objects. In this plot, and in Figures 2-6, the vigorously star-forming galaxy NGC 4449 is 
labeled. Characteristic random errors are indicated separately for the Frei et al. sample and the intermediate- z objects. 

Fig. 2. — Rest-frame U — B versus B — V for the sample samples as in Figure 1. The intermediate-redshift samples of 
paper II largely overlap with the bluest Markarian galaxies, which extend blueward the color-color relation seen for the 
"normal" galaxies from Kent. 

Fig. 3. — Rest-frame i?-band form and scale parameters versus spectral index for the Frei et al. sample. Top panel: 
Average surface-brightness within the half-light radius (SB e ) versus rest-frame B — V. Middle panel: Image concentration 
(C) versus B — V. Bottom panel: 180-degree rotational image asymmetry (A) versus B — V. Characteristic errors are 
given in the top-left corner of each panel. Outlying objects are labeled and discussed in the text. Dashed lines demark 
Early, Intermediate, and Late types in our classification scheme. Symbols are by Hubble type, as defined in the key. 
Different Hubble types are well distinguished, particularly in color. Morphological types are also well separated in C, but 
only the earliest types are well separated in SB e and A. 

Fig. 4. — Rest-frame £?-band parameters of form versus scale for the Frei et al. sample. Top panel: Image asymmetry (A) 
versus average surface-brightness (SB e ). Bottom panel: image concentration (C) versus SB e . Outlying objects are labeled 
and discussed in the text. Dashed lines demark Early, Intermediate, and Late types in our classification scheme. The 
separation of morphological types is less clear than in Figure 3, but the different Hubble types are reasonably segregated. 

Fig. 5. — Form versus form parameters for the Frei et al. sample: Rest-frame _B-band image asymmetry (A) versus image 
concentration (C). Outlying objects are labeled and discussed in the text. Dashed lines demark Early, Intermediate, 
and Late types in our classification scheme. The separation of morphological types is less clear than in Figure 3, but is 
comparable to figure 5 where the different Hubble types are reasonably segregated. 

Fig. 6. — Comparison of form, scale, and spectral index correlations between Frei et al. and Kent samples. Top panel: 
average R band (Kron-Cousins) surface brightness within the half-light radius, SB e (R c ), versus R- or r-band image 
concentration, C(R). Middle panel: SB e (R c ) versus rest-frame B — V. Bottom panel: C(R) versus rest-frame B — V. 
Structural parameters: We measured half-light radius and image concentration for the Frei et al. sample using their R 
or r-band CCD images and circular photometry apertures. Kent measured these structural parameters using elliptical 
apertures on F-band CCD images. Photometric parameters: The Frei et al. sample is subdivided between objects observed 
at (a) Lowell Observatory (filled squares), (b) Palomar Observatory (dotted-circles), and (c) an overlapping subset of the 
Frei et al. sample with existing i? c -band photometry from Buta and Williams (1996; outlined-triangles) . For (a) and (b) 
we used the zeropoints from the Frei et al. image headers (DNATO_BV), and transformations from Thuan-Gunn r and 
Gullixson et al. R to Cousins R c from Frei & Gunn (1994). We have transformed Kent's photometry reported in the 
Thuan-Gunn r-band to R c again based on transformations in Frei & Gunn (1994); Kent corrected surface brightnesses to 
"face-on" values. Regressions: Lines indicate ±lcr about linear least-squares fits to the correlations (dotted, Kent; dashed, 
Frei et al. ) using an iterative clipping method (±2.5cr clip; 10 iterations). In the top and middle panels only the Lowell 
subset of the Frei et al. sample was used in the regressions. The substantial scatter in the Frei et al. SB e (R c ) values we 
infer is due primarily to zeropoint difficulties; we detect no noticeable systematics effects with inclination in SB e (R c ). 
The difference in the correlation between SB e (R c ) and C(R) is largely due to the shallower trend in C(R) with B — V 
for the Frei et al. sample. This may be due to differences between circular versus elliptical apertures. While elliptical 
aperture photometry provides greater dynamic range in C(R), the correlation of C(i?) with B — V has larger scatter. 

Fig. 7. — Form parameters and spectral index for 70 galaxies from the Frei et al. sample as determined by Brinchmann 
et al. Top panel: B band image concentration versus rest-frame B — V. Middle panel: B band asymmetry versus rest-frame 
B — V. Bottom panel: B band asymmetry versus concentration. The asymmetry parameter was determined in a very 
similar manner as our own and thus should have comparable dynamic range. Since our C parameter is logarithmic, we 
plot the logarithm of the Brinchmann et al. C values. These plots are displayed so that they may be directly comparable 
to Figures 3 and 5. The trend in asymmetry for the different Hubble types is more apparent in Figure 3 and 5. In the 
concentration-color plane, the distributions are similar for both studies, although we find a smaller scatter among the 
late-type galaxies, and a larger scatter among the early-type objects. 

Fig. 8. — A representative subset of galaxy images from the Frei et al. catalog, block-averaged by factors 1, 2, 4, and 
6 (top to bottom). While the apparent change in qualitative (visually-assessed) morphology is small, the effects on the 
quantitative parameters C and A can be substantial. Half-light radius and surface-brightness are only weakly affected. 
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Fig. 9. — Resolution dependence of image concentration, C, for the galaxies in the Frei et al. catalog: AC versus the 
half-light radius R e (in pixel units of the block-averaged images). AC is the difference between the concentration index for 
a given simulated value of R e relative to the original concentration value (i.e. that value measured on the observed image). 
Measurements for six definitions of the concentration index are plotted (two types per panel, labeled by line- type). The 
central line (bold) is the median value of this difference, and the bounding lines are the 25% and 75% values, i.e. 50% of 
the simulations are contained between the upper and lower lines for each index. 

Fig. 10. — Wavelength dependence of structural parameters for galaxies in the Frei et al. sample, plotted versus galaxy 
rest-frame color. Dashed lines show the mean differences between blue and red bands and the error bars show the la 
dispersions for three bins in color: B — V < 0.62 (late-type), 0.62 < B — V < 0.87 (intermediate- type), and B — V > 0.87 
(early- type). Top panel: Image concentration Cb — Cr. Nearly all galaxies are more highly concentrated in the red band 
than in the blue, and thus fall below the dotted line at Cb — Cr — 0. This difference is slightly larger for galaxies with 
intermediate- type morphology. Middle panel: Image asymmetry As — Ar. Late- and intermediate- type galaxies are more 
asymmetric in the blue band than in the red band. Red objects are generally very symmetric in both bands, and have 
Ab — Aji ~ 0. This panel can be compared to Figure 2 in Consclicc et al. (1999) where Ab — Ar is plotted versus Ar. 
Since asymmetry and color are strongly correlated for the Frei et al. sample (as seen in Figure 3), the trend in Conselice's 
plot is similar to what is shown here. Bottom panel: Half-light radius i? ej s — Re,R- Although the scatter in this diagram 
is relatively large, it is clear that the half-light radius shows little wavelength dependence over this wavelength range (cf. 
de Jong 1995). Objects of intermediate B — V color tend to be slightly larger in the blue band than in the red, but this 
trend is not seen for either the bluest or the reddest objects. 

Fig. 11. — Axis ratio dependence of image concentration C for galaxies of different morphological types in the Frei et al. 
sample. The dotted line separates the sample into two bins at logio[a/6] = 0.3, corresponding to an inclination of 60°. 
The dashed lines and error bars show the mean and la dispersion for each morphological type and bin. Labeled objects 
are discussed in the text. Top panel: Late-type objects with high inclination have slightly higher measured C than more 
face-on objects. Middle panel: For intermediate-type galaxies the measured concentration indices show no correlation 
with the axial ratio [a/6]. Bottom panel: For early-type galaxies, a tendency can be seen where objects with larger axial 
ratio [a/b] are measured to have lower image concentration. 

Fig. 12. — Axis ratio dependence of half-light radius R e for galaxies of different morphological types in the Frei et al. 
sample. The dotted line separates the sample into two bins at logio[a/6] = 0.3, corresponding to an inclination of 60°. 
The dashed lines and error bars show the mean and la dispersion for each morphological type and bin. Top panel: The 
measured R e are slightly larger for late-type objects with high axial ratios. Middle panel: Intermediate-type objects have 
somewhat smaller R e for high values of [a/b). In both of these panels the scatter is large and the differences between 
the bins are small. Bottom panel: Early-type galaxies with larger axial ratio [a/b] are measured to have <; 30% smaller 
half-light radii. 
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